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ABSTRACT 

We derive formulae connecting the frequency variations in the spectrum of 
solar oscillations to the dynamical quantities that are expected to change over 
the solar activity cycle. This is done for both centroids and the asymmetric part 
of the fine structure (so-called even-a coefficients). We consider the near-surface, 
small-scale magnetic and turbulent velocity fields, as well as horizontal magnetic 
fields buried near the base of the convective zone. For the centroids we also 
discuss the effect of temperature variation. 

We demonstrate that there is a full, one-to-one correspondence between the 
expansion coefficients of the fine structure and those of both the averaged small- 
scale velocity and magnetic fields. Measured changes in the centroid frequencies 
and the even-a's over the rising phase solar cycle may be accounted for by a de- 
crease in the turbulent velocity of order 1%. We show that the mean temperature 
decrease associated with the net decrease in the efficiency of convective transport 
may also significantly contribute to the increase of the centroid frequencies. Al- 
ternatively, the increase may be accounted for by an increase of the small-scale 
magnetic field of order 100 G, if the growing field is predominantly radial. 

We also show that global seismology can be used to detect a field at the level 
of a few times 10^ G, if such a field were present and confined to a thin layer near 
the base of the convective envelope. 

Subject headings: Sun : Helioseismology, solar variability 



1. Introduction 

We study global changes, over the solar 
cycle, in the sun's eigenmode frequencies - 
centroids and asymmetric fine structure - 
in a search of physical changes occurring 
beneath the photosphere. There is abun- 
dant phenomenological information about 
the helioseismic changes, but there is no 
satisfactory physical model describing the 
changes. We consider three possible dy- 
namical sources of the evolution - changes 
in the sub-photospheric small-scale mag- 
netic and velocity fields and a large-scale 
toroidal field buried in a thin layer near 
the base of the convection zone. Here, 
we develop the formalism needed to con- 
nect these dynamical changes to frequen- 
cies changes. 

In our treatment of the small-scale mag- 
netic field, we generalize the method of 
Goldreich et al.(1991, GMWK) to include 
the generalized effect of the small-scale 
magnetic field on nonradial modes, while 
further generalizing to a non-spherical dis- 
tribution of the averaged field. Although it 
is true that radial modes may adequately 
represent lower degree (up to about I = 60) 
nonradial p-modes, if the magnetic field ef- 
fects were confined to the outermost layers, 
this is not true for higher degree p-modes 
or most f-modes. Still, the more impor- 
tant generalization is that we treat a non- 
isotropic, non-spherical field distribution, 
which allows us to interpret the observed 
evolution of the anti-symmetric part of the 
fine structure in the spectrum of solar os- 
cillations (so-called even-a coefficients). 

Furthermore, we study effects a small- 
scale, random velocity field. A role for the 
changing turbulent velocities has been sug- 



gested by Kuhn (1999). However, a first- 
principles treatment still needs to be made. 
We give an estimate of the associated tem- 
perature change and its effect on oscillation 
frequencies. 

Finally, we consider the effect of a 
buried toroidal field, which may be ex- 
pected to be confined near the base of the 
convective envelope. The present work 
represents an advance over earlier ones 
(Gough & Thompson ,1990; Dziembowski 
& Goode, 1991) because wc make a more 
explicit and useful formulation by eliminat- 
ing derivatives of the unknown dynamical 
quantities. This improved development al- 
lows us to obtain more physically revealing 
formulae. The application of this part of 
our work is determining a stringent limit 
on the size of a buried toroidal field. 

2. The Helioseismic Data 

Solar frequency data are usually given 
in the form 

J^Z ^i>in + ^ak,tnVi{m), ( 1 ) 

fc=i 

where the V are orthogonal polynomials 
(see RitzwoUer & Lavely 1991 and Schou et 
al. 1994). The remaining symbols {nim), 
in this equation have their usual mean- 
ings. This representation ensures that P^^„ 
- the centroid frequencies probe 
of the spherical structure, while the a2k - 
the even-a coefficients - are a probe of the 
symmetrical (about the equator) part of 
distortion described by the corresponding 
P2k{cos9) Lcgcndre polynomials. We note 
that in lowest order, perturbations that are 
symmetrical about the equator induce an 
asymmetric change in the fine structure of 
the oscillation spectrum. 
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For the angular integrals, we have 

'■■2tt /•! 



(2) 

where /i — cos 9 and 



Sk,e = (-1) 



, {2k-l)U (2£+l)!! (£-1)! 
A;! (2£ + 2k + (£ - ky: 



Following our earlier works (see e.g. 
Goode & Dziembowski, 2002), we use here 
the following convenient quantities, jk,em 
through the following two relations. 



and 



Ik, in 



(3) 



(4) 



where is the dimensionless mode in- 
ertia calculated for our reference model. 
A clear advantage of the 7's is that their 
growth rephcates the growth of other mea- 
sures of solar activity. For the p-modes, 
the 1// factor takes care of the I- and most 
of the zz-dependence in Az/ and in the even- 
a coefficients. The fact that the residual 
i/-dependence is weak points to a localiza- 
tion of the source of the observed frequency 
changes close to to the photosphere. 

The numerical values of the 7's scale 
with the square of the eigcnfunction nor- 
malization at the photosphere. The nor- 
malization we adopted in our analyses 
of the SOHO MDI data (e.g. Goode & 
Dziembowski, 2002) and which is used 
throughout present paper, is explained in 
the next section. With this normalization, 
the value of 70 reaches up to the 0.3/xHz 
range. The absolute values of 71 and 73 are 
about twice larger. Having determined the 



set of 7fc, one may construct seismic maps 
of the varying sun's activity (Dziembowski 
& Goode, 2002), that is the 7(/u) depen- 
dence. In such maps, a rising 7 reflects the 
local rise of in the activity. The highest 
values of 7(/x) are about 1//Hz and they 
are reached at /i ~ 0.3 and at the peak 
of the activity. At activity minimum the 
highest 7(Ai) ~ 0.2/xHz occurs in the polar 
region. 

In the subsequent sections, we will con- 
nect the 7's, to magnetic and velocity fields 
that are expected to change in the sun over 
its activity cycle. To achieve this, wc start 
from a variational principle for oscillation 
frequencies. In our expressions, the £, n 
subscripts and the m superscript will not 
be given unless it is necessary for clarity. 

3. Variational principle for oscilla- 
tion frequencies 

There are two ways of deriving the vari- 
ational expression for oscillation frequen- 
cies. Both rely on the adiabatic approxi- 
mation, which is adopted throughout our 
study. The first of the two approaches be- 
gins with the linearized equations of fiuid 
motion about a steady configuration (see 
e.g Lynden-Bell & Ostriker, 1967, LEO). 
The other uses Hamilton's principle (see 
e.g. GMWK; Dcwar, 1970). Here, we use 
the form given by LBO with some simpli- 
fication of the variational principle, while 
adding the all-important contribution of 
the magnetic field, as calculated explicitly 
by Dewar (1970). The LBO form is valid 
for strictly steady velocity fields. However, 
we make certain simplifications, which will 
be explained later, to make it applicable 
to statistically steady fields. With this, we 
write 
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u/l = -uC + D 



where 



/ 



(fxp\i\ 



(5) 
(6) 

(7) 



and where v represents the velocity field. 
The eigenvectors, in a spherically-symmetric 
and time independent model of the sun, are 
expressed in the following standard form 

$, = r[y{r)er+z{r)VH]Yr{e, 0) exp(-iu;t). 

(8) 

We adopt some approximations regard- 
ing the eigenfunctions. In addition to adi- 
abaticity, we assume the Cowling approx- 
imation is valid, which is well-justified in 
our application to solar oscillations. Fur- 
ther, we will make use of the fact that the 
oscillations arc cither of high degree or high 
order, which means that 

|£| « Max(r|^;,U|£|). 

Equivalent approximations were also made 
by GMWK but, in addition, they ignored 
the angular dependence of the displace- 
ment. 

Like LBO, we separate the various con- 
tribution to 

D^Dp + Dg + D^ + Dm. (9) 

The pressure term, 

Dp^ J d^xp[E + (r - l)|div£|2], (10) 

is the same as in LBO. The quantity F, 
usually denoted as Fi, is the adiabatic ex- 
ponent. The quantity S is a completely 
contracted double dyadic product. 



: V^. Adopting the standard sum- 
mation convention, we have 

S - Cr,k^k;j = {Cj^k;j - &yO;k + Idiv^p, 

where the subscript ";" denotes covariant 
derivatives. However, with our approxi- 
mation regarding ^, contributions from the 
terms involving the Christoffel symbols are 
negligible, and the derivatives may be re- 
garded as component derivatives. In terms 
of the radial eigenfunctions, y and z, with 
the adopted approximations we have 

S ^ [r(yri/;,-yA);, + A2]|yrl' + 

r{yz),r\VHYn' + [0.5(|/rz^, - zX)\Yn^e 
+z'VHYr-V„Yre];0 + [-h, 



where 



A = y^ - z 



is radial eigenfunction corresponding to 
div^, g is the local gravity, and c is the 
speed of sound. The last term in S is ob- 
tained from the preceding one by the re- 
placement 6 ^ (j). Further, in the adopted 
approximation, we have 

ry.^ = A Az, (11) 



where A = £(£ -|- 1) and 



nVc^ 



rz,r = y-\[- \ -. (12) 
\uj J gr 

The term containing the Brunt- Vaisala fre- 
quency, N, is of the same order as the first 
one for p-modes only, and only in the out- 
ermost layers. However, the whole con- 
tribution from the term involving rz..,. is 
small. Hence, we will ignore the term, so 
that 
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order perturbational expression, which is 



(AiF^r + ivn^rr) + 

i-U + i-h- (13) 

The exphcit expressions for the last two 
terms will not be needed. 

The gravity term simplifies to 

Dg = -2j d'xp^lirW (14) 

after using the Cowling approximation, 
while the velocity term, 

Dv^-J d'xp\{vV)^\^ (15) 

is the same as in LBO, where it was derived 
for a steady field velocity field. We will 
use the same form in our application to 
a statistically steady turbulent field. The 
expression for the magnetic term, which is 
taken from Dewar (1970), is 



D 



M 



1 

47r 



(B-V)^r- 



2div^*B • (B ■ V)C + 

l\B\\E+\diy^\') 



(16) 



We now perturb eq. [5] about the static, 
non-magnetic equilibrium state. The A 
denotes changes in parameters relative to 
this state. However, for centroid frequen- 
cies, A is defined with respect to activity 
minimum because we do not have models 
of the sun predicting centroid frequencies 
with //Hz precision. Since we want to con- 
sider terms that are quadratic in velocity, 
in principle, we need to consider a second 



C 

21 ^ 8P^ 



uAC 



2Iuj 

AD, + D^ + Dm 
2Iu 



(17) 



where 



ADs = A{Dp + Dg) - Lu^AI. (18) 

Actually, we do not calculate the C inte- 
gral or its perturbation, but only comment 
on the role of the terms in C for various ve- 
locity fields. The first term, which is linear 
in velocity, results from rotation, and gives 
rise to odd-a coefficients, which we are not 
treating here. It may be easily shown that 
both meridional and statistically steady 
turbulence do not contribute. The sec- 
ond term, oc due to rotation gives a 
negligible contribution ( Dziembowski and 
Goode, 1992) to p-mode splitting. The 
term oc AC arises from the first order per- 
turbation of the eigenf unctions due to the 
velocity fields. Here the contribution from 
rotation and meridional circulation can be 
shown to be negligible. The only quadratic 
effect of rotation, which we found to be sig- 
niffcant for p-modes is that of the centrifu- 
gal distortion. Thus, it is included in the 
Ds term. For the f-mode even-a's, which 
are not accurately determined, the terms 
involving may be important. The alter- 
native approach, which has been used by 
us in all our analyses of the even-a coeffi- 
cients, is to evaluate the centrifugal contri- 
bution and subtract it from the data. We 
neglect contribution from turbulence to the 
term oc AC, because we include only ef- 
fects of interaction of oscillations with the 
averaged velocity fields. It has to be kept 
in mind, however, that not all effects of 
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turbulence are included in our formalism. 
So that we are left with the expression 



2Iu 



{AD, + D^ + D 



M 



(19) 



In the term, we consider only effects of 
the turbulent pressure and we will be in- 
terested in the part that may vary with the 
solar activity. Wc do not have yet observa- 
tional evidence for changes in turbulent ve- 
locity but such changes are expected. The 
only global changes in velocity which were 
definitely detected are the torsional oscil- 
lations but their effect on frequencies we 
estimated to be insignificant. 

The Ds integral may be calculated con- 
sidering either Eulerian or Lagrangian per- 
turbations. The results must be the same. 
The Dy and Dm integrals are treated as 
perturbations. We may see that the inte- 
grands do not involve differentiation of the 
unknown characteristics of the velocity and 
magnetic fields, rather the differentiation 
is placed upon the eigenfunctions, which 
are known. This is clearly advantageous 
and we will apply the same strategy in the 
evaluations of ADg. 

The calculated frequency perturbation 
for individual (^nm)-modes arc linked to 
the 7's defined in eqs. [3] and [4] by the 
following relation 



with 



27J- ^ 

hn jt=o 

ien = / dxx'^pS, 



(20) 



X — 



p 



p 

— 1 

p 



and 

As a normalization of the eigenfunctions, 
we adopt 

Z/^,n(rphot) = 2 X lOl 

This is an arbitrary choice which leads to 
maximal 7's in the (0.2 — l)/iHz range, 
which is the same order as the frequency 
shifts. This normalization is assumed in all 
expressions for 7 provided in this paper. 

4. Dynamical perturbations of the 
structure 

Here, we include the dynamical effects 
of the magnetic and those of the ve- 
locity fields, V. We write the condition of 
mechanical equilibrium, in the presence of 
perturbing force F, in the following form. 



where 
F = 



Vp + pgcr = . 
dH eg 



dV 
dr 



do r 



1 

47r 



-VB^ - {B ■ V)B 

2 



(21) 

(22) 
p{v ■ V)v 



We neglected the perturbation of the 
gravitational potential, which is justified 
because we considering perturbing forces 
concentrated in thin layers containing lit- 
tle mass. eq. [21] implies 



and 



p = -7i + h{r) (23) 

The quantities V and Ti. represent non-gas 
pressures, which in general are anisotropic. 
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It should be noted that when the non-gas 

pressure is isotropic, then the mass dis- 
tribution remains spherically symmetric. 
The quantity h{r) may only be determined 
by utilizing the condition of thermal equi- 
librium. For the non-spherically symmet- 
ric parts of the force, the pressure and den- 
sity follow from the condition of mechani- 
cal equilibrium. 

From now on, we treat F small per- 
turbing force. Primed letters denote Eule- 
rian perturbations of the respective struc- 
ture parameters, letters preceded by 6 de- 
note Lagrangian perturbations, and letters 
without such symbols imply unperturbed 
variables. We use the standard relation 



and from mass conservation 



■5/ 



dr 



and we adopt SMj. = 0, for both spherical 
and aspherical perturbations. 

Note that if we make a Legendre expan- 
sion in even orders, P2k{cos9), of H and 
V, all the expansion coefficients and p'^, 
starting from k — 1 are completely speci- 
fied. Here we are considering only even or- 
der polynomials because these are the ones 
that contribute to the even-a coefficients. 
The anti-symmetric (odd-order) polynomi- 
als average out. From eqs. [23] and [24], we 
get for k > 0, 



P'k 



and 



p'k 



1 d 

g dr 



—'Hk 



{Hk — Vk 



(25) 



The expansion coefficients for the La- 
grangian perturbations are calculated as 
follows. From the radial component of 
eq. [21] we have 



SkP = -Vfe + 4 



/ 



-gpdr fti -Vk (26) 



dr 



SkP _ rt ^kr 



SkP 



(27) 



The approximate equality in the pre- 
ceding equation corresponds to neglecting 
of the perturbation in the mass distribu- 
tion above the point under consideration. 
This is certainly valid for all of our applica- 
tions. The approximate equality in eq. [27] 
is just the local plane parallel approxima- 
tion. This is valid for most of applications 
considered here. The only possible excep- 
tion will be discussed briefly in subsection 
5.1. Both approximations were adopted in 
GMWK. We stress, however, they are not 
needed for deriving expressions for 7^, ex- 
cept for k = 0. For k > 0, it is only impor- 
tant and, in fact well justifled, for seismic 
determination of the aspherical part of the 
subphotospheric temperature changes. To 
this aim, we first derive an expression for 
Skr from the relation between p'^ and SkP, 



SkT = 



9P 



(28) 



Then, using the linearized p{p, T) relation 
we obtain 



s_kr 

T 



XtP 



Xp 



d 



d\np 



+ 



(29) 



Here, we used a standard notation in astro- 
physics, e.g. x's denote derivatives of logp 
with respect of the log p and log T. We see 
that for the non-spherical part, all pertur- 
bations of thermodynamical quantities are 
determined by H and V. This is not true 
for A; = 0, where one of the thermodynam- 
ical parameters is left free. Choosing the 
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temperature, we have 

= — + XT 

P Xp \P 



T 



(30) 



or, if we choose the entropy per mass vari- 
ation, ^o^*, in place of ^o^; 



^oP 1 Vo XT SqS 



(31) 



p r p Xp cp 

4.1. Turbulent pressure 

The large-scale average of the Reynold's 
stress, F — —p{v ■ V)^, due to the tur- 
bulent velocity, is evaluated in the local 

Cartesian system with axes parallel to e = 
(cr, eg, e^) (effects of curvature are neg- 
ligible for this small-scale velocity field). 
Then, we have 



Fi = -{pVjVi).. + {pVj).jVi. 
We use the following relations 

and 

dp 



{pvj).jv, = -—Vi = 0, 

where double-subscripted S is, as usual, the 
Kronecker symbol. That is, we assume 
uncorrelated velocity components and the 
rate of density fluctuations. Hence, we 
have 



Fi — —{pvf).- no summation over i ! 

(32) 

Further, we allow the vertical component 
(r) to be statistically different from the two 
horizontal {6 and 4>) ones. In this, values of 
pVj are treated as functions of depth, and 
slowly varying functions of the co-latitude. 



The latter dependence is represented in the 
form of a Legendre polynomial series. 



pviVj = pSijY^[SjrT,y {r) + 
1 



k=0 



-T,"{r){6,e + 6,^)]P,,{cose), (33) 

where wc included only terms that are 
symmetric about equator. 

Inserting this expression into eq. [32] 
and using the definition of V and Ti as 
given in eq. [21], we get 

Vk = pT^, and Hk = ^pT,^. 

(34) 

4.2. Small-scale random magnetic 
field 

Our treatment of the small-scale mag- 
netic field is analogous to that of the tur- 
bulent velocities. That is, the correlation 
matrix for the field components is repre- 
sented in the form of the following Legen- 
dre polynomial series. 



1 



fe=0 



-M^(r)(5jg + 5j^)]P2k(cose). (35) 

Components of the Lorentz force, treated 
locally as Cartesian, are given by 

Averaging over wide zonal areas and mak- 
ing use of Bi-i = 0, we get 



BfjBr-g + B^Br-, 



IdBl 



8 



Thus, 



J_d_ 
Stt dr 



which is the same expression that was ob- 
tained by GMWK. Let us note that the net 
effect of the vertical component of the ran- 
dom field magnetic on the vertical struc- 
ture is opposite to that of the horizontal 
components. The radial component acts 
as a negative pressure, when it rises the 
gas pressure must rise too. 

To evaluate the horizontal force, we use 



to get 



B^Bn-d, — -{Bg).0 + Br-rBe, 



1 

4:71 L 



{Br Be)., + 



Bl - Bl 



Stt 



Finally, for the coefficients in the expan- 
sion of the magnetic pressure, we obtain 



Ml 



V 



Stt 



and Hk — 



Ml 
Stt 



(36) 

Just as in the case of the turbulent ve- 
locity field (see eq. [34]), isotropy implies 
Vfc = Ti-k, hence no density perturbation, 
but effects of a departure from isotropy are 
clearly different. 

4.3. Large-Scale Toroidal magnetic 
field 

The large-scale field, B = Bt{r,9)e^, 
gives rise to the Lorentz force 



F = 



V {By sin'' 9) V{Bfsm''9) 



Snr'^ sin^ 9 



STTsin^^ 



The last approximation is valid for a field 

confined to a narrow layer, which we will 
assume here, so that we have 

Stt 



and 



d9 d{B^sm'9) 



sm^9 



d9 



We now put Bt (r, 9) in the form of the 
following series. 



Bt{r,9) = J2 



2j + l ^ . 

27(JTi)^*'^-^^- 



dPj{cos9) 
d9 



(37) 

Note that with this representation, Btj is 
the surface averaged intensity of the field 
component at a distance r from the cen- 
ter. Considering only first two terms in the 
expansion, we get the following non-zero 
components of the Legendre polynomials 
expansion for V 



Vi 



and for 7i 



1 



Bt,i + -Bt,2 



IGtt 



(38) 



167r 

V2 



B, 



\i + ^Bi^y 



" 28n ' 



(39) 



SGtt 



(40) 
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5. The term arising from the struc- 
tural perturbation, ADg 

The frequency perturbation arising through 
the perturbation of the structure for all 
forces considered by us is, typically, of the 
same order as that arising directly from 
the forces. We will express now perturba- 
tion the structural term Dg in terms of V 
and H calculated in the previous section. 
The variational principle ensures that we 
may keep ^ (not y and z !) unperturbed. 

5.1. Calculations of ADg for cen- 
troid frequency shift 

Here, using the Lagrangian formulation 
of the perturbations is more convenient. 
Since we have S{drr'^p) = (we drop from 
here on the subscript at S if it is zero) , there 
is no contribution from A7. Furthermore, 
with our approximation for the eigenfunc- 
tions, the contribution from Dg is negligi- 
ble. For the present application, it is con- 
venient to write eq. [10] in the form. 



Dp — J drr^p ^ 



dZ 
dr 



where Z = 2r{Ayz — y"^) and, which after 
integration by parts becomes 



drr 



gp-^]z+pry 



In the whole solar envelope the second 
term in the coefficient at Z is much less 
than the first one and it will be ignored. 
Now we calculate A.Ds,o ~ A.Dp using 



Si 
9 



SZ 



-2- 



5r 



and 



^(i + r,)-^(i-r,). 



where we denoted by Fj, and Fp logarithmic 
derivatives of F. Further, we use eq. [26] to 
eliminate 5p and eq. [30] to eliminate 5p. 
Finally, with our approximations regarding 
the eigenf unctions, we get 



where 

^isoth — 



I 



p 



drr 



AsothVn + 



^ ST ^ 6r\ 



(41) 



F 

2Az\ 



^ Xp > 



+ 



X/0 ^ 



(42) 



XT, 



= -^F[(l + Fp)A2 + 2A^A], (43) 
Xp 



and 

Vr = 2pFA^A + 6— (t/2 - Ayz). 
P 

The relative roles of the temperature 
and radius depends on the character of 
perturbation and mode. As pointed out 
by Dziembowski, Goode and Schou (2001, 
heretoforward DGS), the latter may be- 
come dominant for f-modes, if the mag- 
netic perturbation is predominantly below 
the region sampled by these modes. For 
f-modes, to a very good accuracy, we may 
use 

(44) 



V, = -Q^Ayz 



P 

In Section 9 of the present paper, we will 
discuss in greater details the role of tem- 
perature and radius variation in the f- and 
p-mode frequency changes. 

If instead of eq. [30], we use eq. [31], then 
we get an alternative expression for AD^ o 
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which is particularly useful if the perturb- 
ing force is localized in the deeper layers, 
which may be regarded adiabatic on the 
eleven-year scale, 



AD, 



drr 



, ^ 5S ^ 5r\ 
p ( Vt— + Vr— 

Cp r J i 



,(45) 



where 

Pad = -[r(i+rp)+i+r^]A'-2A^A. (46) 

5.2. Calculations of ADs for the 
splittings 

In the present application, it is more 
convenient to treat the perturbations of 
the structural parameters as being Eule- 
rian. We consider distortions proportional 
to P2k- We will see that within our approx- 
imation, all the angular integrals appear- 
ing in A/, ADp and ADg reduce to Qk- 
These factors take care of the k and m de- 
pendence. The property is self-evident in 
the case of Dg. From eq. [14] with the use 
of the definitions given in eqs. [8] and [2], 
we get 

ADg = -2Qk j drr^gp'^y'' 

and with eq. [25] after one integration by 
parts, and use of eq. [11], we get 

ADg = AQk j drr\Vk - nk)y{X + Az). 

(47) 

The cases of Ip and Dp are more involved. 
We first note that 



I d{cos9) J 



AQk. 



The approximation assumes i ^ k, which 
is not valid for low degree modes. However 
the terms involving this factor are signifi- 
cant for such modes only in the core, which 
we assume is unperturbed. Thus, for A7 
we have approximately 

AI = Qk J drr^p'^S. 

Again, we make use of eq. [25] and inte- 
grate by parts, and with eqs. [11] and [12], 
to obtain approximately 



AI = 2Qu 



dr-{Vk 
9 



Hu){y\ + 2Azy). 

(48) 

In calculating ADp, we first note that the 
term in eq. [13] does not contribute, 
which follows from the assumed axial sym- 
metry of the perturbation. The contribu- 
tion from the {--^e term is nonzero, but 
it is small, as may be justified as follows. 
Integrating by parts over 6, one gets the 
k{k + l)Qk factor from the angular integral 
and the whole contribution from this term 
is of the same order as the one neglected 
above. Thus, we have 



ADr. 



Qk j drr' [r(p;(i + r^) + 



The quantity is again eliminated by 
integration by parts. The use is also made 
of eqs. [11] and [12]. The result is 

ADp = -Q, J drr^nk[Til + Tp)\'' + 

2A{s + zX)] + {nk-vM, 



where 
^ = Fp 



( ^ (i In \ , n 



o?lnp 



11 



2\{Kz - —y 



dlnp 



Combining this with eq. [48] in eq. [18] 
(ADg is neghgible), we obtain 



AL>,,fc = gfc J drr^iV^k+VfUk), (49) 

(50) 



where we denoted 



c 



2 

uj r 



9 



{yX + 2Azy), 



and 



= 2C - V - r(l + rp)A2 - 2A{£ + zX). 

(51) 

6. Frequency change to due varying 
turbulent pressure 

For evaluating Dy according to eq. [15], 
we use the random velocity field represen- 
tation given in eq. [33]. We note that 

k 

lr,"A")p,„ (52) 



where 
and 



dy 
dr 



m\2 



+ r 



dz 
dr 



IVY, 



m|2 



A 



H 



y^\vYr\'' + z^wixn-fi-fi? 

The radial derivatives in are eliminated 
with the help of cqs. [11] and [12]. The an- 
gular integrals are evaluated by parts keep- 
ing only derivatives of the spherical har- 
monics. This approximation justifies, in 
particular, the replacement 



In this way, we get the contribution to 
from the P2k component of the turbulent 
pressure. 



D 



v,k 



Qk J drr^p r/(A^ + 2A^A + 

(53) 



^ 2 



The contribution to Dg^ from the induced 
change in the gas pressure at constant tem- 
perature and radius is given by 

AD, = J drr''pVi,^ti,T^, 

which follows from eqs. [34] and [41]. Using 
these two expressions in eq. [19], we get 



where 



1 



isoth 



2ujI 



(54) 



Asoth - (A' + 2A^A + AS) 



and 



The expression for Asoth is given in eq. [42] . 
The complete expressions for the Tfs 
are given in the Appendix (eqs. [Al] and 
[A2]). Here, we provide only the asymp- 
totic forms of the TZ^s for p-modes where 
i\z/y\ -C 1, which is valid sufficiently 
above the lower turning point, as well 
as, being the form appropriate for the f- 
modes. Our approximation for p-modes 
is the same as that used by GMWK, and 
that made for f-modes is the same as made 
by DCS. 

For the p-modes, the leading terms are 
those proportional to A^. If we keep only 
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these terms, and, in addition, if we ignore 
the derivatives of F, we find 



i;,isoth 



r 



(55) 

For the f-modes. wc have (e.g. Appendix 
in DGS) y £z and \ <ti £z which imphes 



7?V _ 

'%,isoth ~ 



(56) 



Thus, in both hmiting cases a decrease in 
the turbulent pressure results in a frequency 
increase. In fact, this property is vahd for 
all observed solar oscillations. With the 
help of eqs. [45] and [46] , we may easily 
obtain expressions for the adiabatic (ad) 
instead of the isothermal kernels. 

For A; > 0, we get from eqs. [49] and [34] 



1 



This together with the eq. 
eq. [19] gives 



used in 



(Aa;)„,jfc 
where 



Qk 

2hj 



(57) 



and 



V 



(A^ + 2AzX + AS) 



-iV^-AE). 



The general expressions for and 
are given in eqs. [50] and [51], respectively. 
The complete expressions for the 7^'s are 
in the Appendix (eqs. [A5] and [A6]). Note 
that 1^^. and TZ^^. are the same for all k > 
0. 



The asymptotic expressions for p-modes 



are 



7?^ ~ 



where 



-2c-A^ 



c 



'^v,k ~ 



C-^A^ (58) 



2 

00 r 



9 



-yX 



is the highest order term in the cj — ^ cxd 
asymptotics. However, for solar p-modes 
in the outer evanescent zone, C is compa- 
rable to A^, and below it changes from -|- 
to - and therefore we keep terms involving 
both quantities. 

For the f-modes, we now have 



'<'v,k ~ 



-SA£, 



n. 



H 

v,k 



-AS. (59) 



Once wc have the TZ kernels, we can evalu- 
ate the 7's introduced in eqs. [3] and [4] for 
a specified turbulent velocity field. We give 
here an expression, which is convenient in 
application to solar data 



^v,k 



d 



"^phot 



1 Mm 



^v,k 



H 



^''^ V 1 km^s-^ 



1 km^s ^ 
// Hz, (60) 



where (iphot is the depth beneath the pho- 
tosphere, 

<f = 3.76 X x^-'^ljf^^-'pv.lf . 

The normalization of the cigenfunctions in 
7?.'s must be the same as in / used in the 
definition of 7. In the Appendix, we give 
exact expressions for TZ^'^^ . 

In Fig.l, we show examples of the 
kernels that arc important for evaluating 
the p-mode 7's due to the perturbation of 
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0.5 1 1.5 




Fig. 1. — Kernels for calculating 7^^^ ac- 
cording to eq. [60] at three selected fre- 
quencies plotted as functions of the depth 
in outer part of the standard solar models. 
The three kernels are ^-independent for p- 
modes in this part of the sun. The K^q 
kernel is oc A and is much smaller. 



the turbulent pressure. Although the gen- 
eral trend is consistent with the asymptotic 
formulae (eqs. [55] and [58]), there are vis- 
ible small-scale structures arising from the 
derivatives of F, which we have ignored in 
these two asymptotic formulae. The ker- 
nels for multiplying the vertical compo- 
nent, have significantly larger absolute val- 
ues and are negative. Thus, we expect a 
rise of the 7's with a decrease of turbulent 
velocities. Since the increasing magnetic 
activity is expected to inhibit turbulence, 
the trend of the calculated effect in 7 is 
consistent with observations. 

With the help of Figure 1, we may 
roughly estimate the required change in 
the mean turbulent velocity needed to 
account for the measured 7's, under the 
assumption that this change is the only 
source of the 7's. From numerical simula- 
tions, we know (e.g. Abbett et al., 1997) 
that velocity fluctuations at the level of 1 
km/s persist over the whole layer shown in 
figure, with a maximum of nearly 3 km/ s at 
c^phot ~ 0.1. The frequency averaged value 
0/70 ~ O.S/iifz (DGS) requires a fraction 
(0.2 - 0.5) of one percent decrease in the 
radial component of velocity fluctuations. 
The largest (k = ^ and 3) require about 
one percent decrease. Such small changes 
would not be easy to detect. 

Fig. 2 shows that the kernels for the 
f-modes are very different from those for 
the p-modes. In this case, the asymptotic 
expressions for TZ are quite accurate. The 
kernels scale as A/z/ oc t^'^^ . All the kernels 
have similar shapes, and they all negative. 
The value of the f-mode kernels are compa- 
rable to those of the p-modes. The maxi- 
mum measured values of 70 for f-modes are 
about twice that for p-modes, but the er- 
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are transformed as follows 
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Fig. 2. — Kernels for calculating 7^ ^ for 
f-modes at three selected ^-values. 

rors are large. The uncertainty for /c > 
is even higher. 

More detailed analyses of the 7jk(i^) de- 
pendence reflecting the kernels frequency 
dependence seen in Figs. 1 and 2 are 
needed to say more about the nature of the 
required change in the turbulent velocity. 
However, at this stage we may already con- 
clude that this change must be regarded as 
important, perhaps the dominant contrib- 
utor to solar p-mode and f-mode frequency 
changes over the activity cycle. The sign 
of the observed changes agrees with the ex- 
pected inhibiting effect of the field on con- 
vection. 

7. Frequency change due to varying 
small-scale, necir-surface mag- 
netic field 

With the random field being described 
by a single P2fe-component (see eq. [35]), 
the three terms in the integrand of eq. [16] 



M^iX^ + 2AzX + AS) + Mj^^S 



Qki 



-2div|*B ■ {B ■ V)i 
-[2M]^X{X + A^) - M^AzX]Qk, 



and 



Idiv^p 



2' ' ^ 

{Ml + M^)[X' + A{zX + £)]Qk. 



The transformations use integration by 
parts over the surface, our approximations 
regarding the eigenfunctions, and the an- 
gular dependence of the averaged fields. 
Note that the first term is fully analogous 
to the integrand in Dy that was considered 
in the previous section. 

With the above expression, we get from 
eq. [16] 

DM,k = ^ / drr^{SMlA{zX + 28) + 
8M^[X^ + A{2zX + l.h8)]]. 

For the spherically symmetric part of 

ADs, we use eq. [41], ignoring here again 
the temperature and radius changes. This 
combined with eq. [36], yields 

AAsoth = ^ y drr'V^iM^ - M/). 

(61) 

Using last two expressions in eq. [19], we 
get 



M,isoth — 



167ra;/ 

iV 



+^M,isoth-^o'^)) 



M,lsoth-^o' 



(62) 
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where 



'^M,isoth 



2K{zX + 28) - Pi,oth 



and 



'^M,isoth 



2A' + k{AzX + 3£) + Asoth, 



with I>s,isoth is given in eq. [42]. The com- 
plete expressions for the 7^'s are in the Ap- 
pendix (eqs. [A3] and [A4]). The asymp- 
totic expressions for the p-modes are 



7^ 



V 

M,isoth 



r + 



and 



'^M,isoth ~ 



2-T-—)x\ 
Xp) 



(63) 



The adiabatic kernels, equivalent to those 
found by GMWK, are obtained by replac- 
ing Xp with r~^. In both cases the equa- 
tions imply that an increase in the verti- 
cal component leads to a decrease in the 
frequencies, while the opposite is true for 
the horizontal component. The sign of fre- 
quency shift due to the horizontal field is 
opposite to what one might have naively 
expected because the dominant effect of 
such field arises through the perturbation 
of the equilibrium structure (the Ds term) 
and not by the direct effect of the field 
on oscillations (the Dm term). The for- 
mer term is negative because the horizon- 
tal field causes a local expansion, hence 
an increase of the sound propagation time. 
The vertical field has an opposite effect. 
An isotropic {Mq — 2M^) field increase 
implies a net frequency increase but the re- 
quired increase to account for the observed 
frequency changes is large. 



For the f-modes, the X>s,isoth may be ne- 
glected and we have 



'^M,isoth ~ 



M,isoth 



3 



(64) 



Thus, an increase of either component of 
the magnetic field implies a frequency in- 
crease for f-modes. 

For A; > we have from eqs. [49] and [36] 

^ j drr\V^ - VX)Ml 



AD 



and 



J drr\n\,Ml 



(65) 



where 



and 



H 
s 



+ 2A{zX + 28) 



^M,fe-^; +2A^ + A(4^A + 3£:). 

Expressions for T>Y and T>f are given in 
eqs. [50] and [51], respectively. The com- 
plete expressions for the 7^'s are in the Ap- 
pendix (eqs. [A7] and [A8]). Again kernels 
T^M,k are the same for all /c > 0. The The 
asymptotic expressions for p-modes are 



'^M,k ~ 



4c - ^A^ 



7^ 



H 

M,k 



-2C + 2A^ 
(66) 



For the f-modes, we now have 



n 



V 

M,k 



6A8, 



H 

M,k 



A8. 



(67) 



Now we have for 7's 

C^phot 



lM,k = 



d 



1 Mm 



K 



M,k 



+K 



M,k 



5Ml 
1 kG 



1 kG 

H Hz, (68) 



16 



where 



In Figs. 3 and 4, we show kernels for cal- 




0.5 1 1.5 



Fig. 3. — Kernels for calculating centroid 
frequency shifts due to a small-scale mag- 
netic for p-modes at the three selected fre- 
quencies. 

culating jm according to eq. [68] . Note the 
strong sensitivity to the frequency, which 
emphasizes the probing power of the 7(2/) 
dependence. Further, note that the kernels 
imply that an increase in the radial field in 
outer layers will lead to an increase in the 
mean frequency, while that of the horizon- 
tal field has the opposite effect. The growth 
of the vertical field also leads to an increase 
of 7's at A; > 0, but the effect of the hori- 
zontal field growth is impossible to predict 
as it depends a lot on the depth where it 
takes place. Also in the case of magnetic 
fields the kernels for f-modes differ signif- 
icantly from those for p-modes as we may 
see comparing Fig. 5 with those in Figs. 3 
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Fig. 4. — The same as Fig. 3 but for cal- 
culating the even-a part of the frequency 
splitting 

and 4. For f-modes, the effect of the hori- 
zontal components of the field is similar to 
the vertical ones but smaller. 

Again, we may use the plots shown in 
these figures to assess the required mag- 
netic field changes needed to account for 
the measured 7's. Let us first consider 70. 
For p-modes, the minimum requirement 
for the field increase is obtained, if we as- 
sume that only radial component increases 
and it is ~ 100 G (DGS). The number rises 
to above 200 G if we assume an isotropic 
field increase (GMWK, DGS). This latter 
value is unacceptably high. Also, there 
is a higher requirement to account for f- 
mode 7o's. Though the observational ac- 
curacy is poorer than in case of p-modes, 
this may be regarded as a piece of evidence 
against the direct effect of a changing mag- 
netic field as the sole cause of the frequency 
changes. Also accounting for the even-a co- 
efficients sets more stringent requirements 
on the near-surface magnetic field, which 
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Fig. 5. — Kernels for calculating the cen- 
troid shifts (upper panel) and splittings 
due to vertical components of the field for 
f-mode modes at three selected degrees. 
The kernels due to the horizontal compo- 
nents have similar shapes but smaller (fac- 
tor ft! 3 at A; = and ft! 6 at /c > 0). 

may be difficult to reconcile with the mea- 
surements. 

8. Temperature and radius varia- 
tion 

GMWK were first to considered the role 
of temperature variations in the p-mode 
frequency changes. They correctly ob- 
served that the temperature decrease at 
constant pressure results in frequency de- 
crease because the effect of local expansion 
exceeds that of sound speed increase. How- 
ever, they excluded the effect of tempera- 
ture change as a primary source of the mea- 
sured frequency changes. Here we recon- 
sider the effect using our formulation pre- 
sented in Section 5.1. With eqs. [20],[19], 
[41], and [43] we obtain the following ex- 



pression for the temperature contribution 
to 7o. 



7t,o = 



d 

ST 



^phot \ 

1 Mmy/ 

yUHz, 



Kn 



(69) 



where 



and 



-0.72- 



1 mHz \ujr/ 



— I X p — IZi 



Xp 



= (1 + rp)A2 + 2AzX. 

For p-modes the first term is dominant in 
TZt- Except for our taking into account 
the derivative of F, it is the same as in 
GMWK. For f-modes the second term is 
much greater but the entire kernels are 
much smaller than for p-modes, as we may 
see in Fig. 6. Thus, we will consider the 
effect of temperature only for the p-modes. 
With the plots in the upper panel, we may 
estimate that the fractional temperature 
increase in the outer layers at a 10"^ level 
implies decrease of 70 at a 10~^/iHz level, 
which is significant. The question arises 
whether such temperature changes during 
the solar cycle are feasible. 

Temperature variation correlated with 
magnetic variations are expected. How- 
ever even the sign of it is a matter of de- 
bate. Gray & Livingston (1997) put for- 
ward evidence that there an increase of T^s 
between the activity minimum and maxi- 
mum by some 1.5K, that is STes/Tes ~ 
2.6 X 10~^, which would account for the 
observed variation in the solar constant. 
Since the optical depth increases with tem- 
perature increases, 5T at dphot = is some- 
what greater than ^Teflf- An estimate us- 
ing the Eddington approximation yields 
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Fig. 6. — Kernels for calculating the fre- 
quency shifts due to temperature increase 
for p-modes at the selected frequencies 
(upper panel) and f-mode modes at the se- 
lected degrees (lower panel). 

8T{d^i,ot) = 1.45reff. The result of Gray 
& Livingston is not generally accepted. 
Spruit (1991) argues that the dominant ef- 
fect of the magnetic field on temperature is 
through inhibition of convection and hence 
it implies cooler layer outer layers at high 
activity. If this indeed the case, then the 
induced temperature variation contribute 
to frequency increase. Spruit (1991) ex- 
plains the irradiance increase correlated 
with the activity as a result of an increased 
corregation of the photosphere. 

In a crude manner, the expected tem- 
perature change may be linked to the 
change in the turbulent velocity. In Sec- 
tion 6, we found that (2 — 5) x 10~^ change 
in turbulent velocity suffices to explain 
the maximum value of 70 ~ 0.3;uHz. Our 
aim is to estimate the values of 5T /T in 
the subphotospheric layer extending down 



to (say) 1 Mm associated with such a 
change in the velocity. To this aim, we 
rely on the mixing-length approximation 
(MLT) and we mimic the inhibiting ef- 
fect of the field by varying the MLT pa- 
rameter a. While perturbing a, we keep 
both R and entropy constant in the adia- 
batic part of the convective zone. Adopt- 
ing 5v/v = -3 X 10-3, we find 5T/T = 
-1 X 10-^ -2 X 10-^ and - 1 X 10"^ at 
c^phot = 0, 0.5, and 1 Mm, respectively. 
The implied 4 K decrease of Tphot between 
solar minimum and maximum seems un- 
acceptably large. This, however, should 
not be regarded as a case against changes 
in turbulent velocities being the primary 
source of the frequency changes because 
our treatment of energy transport was very 
crude indeed. Rather, we want to empha- 
size here that temperature changes in the 
subphotospheric layers must he considered 
as significant contributor to the observed 
frequency changes over the solar cycle. 

The aspherical part of the temperature 
perturbation is fixed by the condition of 
mechanical equilibrium. Eq. [29] expresses 
5kT in terms of the expansion coefficients 
Vfc and which in turn are linked to the 
expansion coefficients for turbulent pres- 
sure (eq. [34]) and magnetic field (eq. [36]). 
We may see that any inference on temper- 
ature depends on the derivative of the per- 
turbing force, and thus reqiiires a detailed 
analysis of the 7(2/) dependence. Currently 
available data are probably not accurate 
enough for this. 

In contrast to the effect temperature, 
which may only be important for p-modes, 
the effect of radius perturbation is likely 
to play a role only in f-mode frequency 
changes. Considering in eq. [41] only the 
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effect of radius perturbation and adopting 

the approximation 2hyz C,£^ which is 
valid for f-modes, we get from eq. [19] 

This expression was used by DGS, who 
argued that the part of the frequency in- 
crease which is proportional to v may be 
explained by the part of 5r which is com- 
mon to all modes in the f-mode set. Their 
set contained modes with £'s from 137 to 
300. The common part must originate 
below the outer part of the sun sampled 
by all these modes, that is below radius 
r = O.988-R0. They argued that its likely 
cause is an increase in the radial compo- 
nent of the magnetic field by few kG. It 
is ironic that the best evidence for deep 
seated magnetic field changes may come 
from modes which do not directly probe 
the region where the field is located. Un- 
fortunately, what we get with these modes 
is only an integral constraint on the field. 
Of course, it would be advantageous to 
have a direct probe for the deep seated 
field. 

9. Frequency perturbation due the 
horizontal field in deep layers 

It has been argued (see, e.g., D'Silva 
& Howard 1993) that a horizontal field of 
i? ~ 10^ G is present in the region near 
the base of the convective envelope. Seis- 
mic evidence for the presence of such a field 
is still controversial. 

First, we consider a large-scale toroidal 
field of the form given in eq. [37] , and trun- 
cated at j = 2. The consecutive terms in 
L>M(see eq. [16]) are calculated under the 



same approximation as used in §3. The 
three integrands become 

-2diYCB ■ {B ■ V)| = -2Wm^\YpfzX, 
and 

i|s|2(s+|div^n = w{i-f,')\Yr\' 

[X^ + A{S + zX)], 

where we denoted 

Note that the last term in W cancels out 
upon integration. Calculating the surface 
integrals first two terms in Dm-, we rely on 
the following recursion relation (DG 91). 

2h — \ 

= M4A + l-4F) ^^^-^^'^~ 
{2k - If - 2m' + qk-2{'2k - 3){k - 1)}, 

where 

qk = j\^ij\<l)fi^''\Yn^. 

With this relation, assuming A ^ 1 and 
using explicit expressions for the sur- 
face integral becomes 

/I /'27r 
dji J d(t)n'^\Yp\^ ^ 




We could assume A » 1 because for low 
degree modes the third terms dominates. 
The surface integral in this term is easily 
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expressible in terms of Qi and Q^- Com- 
bining all the three terms in eq. [16], we 
obtain 



D 



M 



1 

47r 



drr'{Bl[^{2X' + 3AS) 



~{X^ + 3AS-3AzX) 



^(2A2 + 3A£:) + ^(A^ + 3AzA) 

2g 



^(3A2 + 7A£-5A^A)]} 

We now proceed to calculate the contri- 
bution to the ccntroid frequency changes 
due to the induced adiabatic pressure 
change. The adiabatic approximation is 
now well-justified on the grounds that the 
layer where the field is expected is located 
deep enough. Setting t^S" = in eq. [45] 
and using eq. [38], to express Vq, we obtain 



s,0 



IGtt 

Eq. [46] gives an explicit form of I^s,ad in 
terms of the eigenf unctions. 

The corresponding contribution to the 
splittings is obtained from inserting eqs. [39] 
and [40] into eq. [49] yields 



1 



-^(2Pr+3Pf)4.2 

Combining all three D integrals into a sin- 
gle expression for the frequency shift due 
the j*'^-component of the toroidal field, we 
have 

1 i 



Stt/u; 



A;=0 



drr^Bljnt,ky 



(70) 



The j^^ component generates all the 7's 
from A; = up to J. From the D's cal- 
culated above, we get for the following ex- 
pressions for the 7^'s at j = 1 and 2 

nt,oi = 7^t,o2 = + A" + ^K£, 



1, 



7^t,l2 = - 2Pf + 2 A' + 6AzA) 



7^ 



t,22 



4A(7£ - 5^A)] 



The quantities I>ad, , and , are 
given in eqs. [46], [50] and [51], respec- 
tively. In these equations, one may make 
use of yet another approximation that is 
valid for p-modes in the acoustic propaga- 
tion zone and beneath. For solar p-modes 
(the f-modes are irrelevant here), we have 
A ~ —Ar]z, where 



V = 



1 /ujr^^ 
A \~) 



is the square of the ratio of the mode to 
Lamb frequencies. At the inner turning 
point, we have rj — 1. With the above 
expression for A, we will obtain a more ex- 
plicit form of the TVs. The products of the 
eigenfunctions occurring in the 7^'s may 
now be expressed as follows 

zX = -7^£^, A^ = Ary2£^, £^ = Az\ 

(71) 

C = y^(2 - ^) = —A'^r]{2 - r])yz. 

9 gr 

(72) 
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Finally, we assume the ideal gas law, which 
is a fully adequate approximation in the 
region considered, to obtain 



7e, 



tm 



i,02 



A 



12 V 3 ' 



7^t,ll ~ -C - A 



t,12 



15^ 5^ 



2S + 7](^^7i-5^ 



7^t,22 ~ -y - ^[lo^ - ^(3r? + 38)£:^] 

We now rewrite eq. [70] in the following 
convenient form, specialized for the sun 



7,. = e/ ^-^^-'^^ (nib) 



j>k 

with X = t/Rq 



//Hz, 
(73) 



K,,,, = 1.48 X 10-^l^x'nt,k. 



The modes that are most sensitive to 
the field in the vicinity of the base of 
the convective envelope are those of mod- 
erate degree with turning points located 
there. This is illustrated in Figs. 7 and 
8, where we show the kernels Kt,k{x) for 
three £ = 30 modes. The modes have fre- 
quencies in the 1.95 to 2.63 mHz range. 
The lower turning points range, accord- 
ingly, from X = 0.742 to 0.654. The n=6 
mode, which has its inner turning point at 
X — 0.721 that is above the base of the 
convective zone, probes the field not only 
within the convective envelope, but also 
the region immediately beneath, which is 



below its inner turning point! This lat- 
ter fact is in contrast to the ray approx- 
imation in which this mode would know 
nothing about the region beneath its inner 
turning point. The turning of the n = 8 
mode is at x = 0.68 and this mode is the 
best probe of the bulk of the overshoot- 
ing zone that extends down to about 0.70. 
Similar results would apply for different i- 
values, after an appropriate re-scaling of 
the frequencies so that the r] (x i/u ratio 
is preserved. We see that the toroidal field 
increase leads to a corresponding increase 
in 7o and that the effect is mostly of the 
opposite sign for A; > 0. 

If a toroidal field of 1 MG would pre- 
vail in the layer between x — 0.68 and 
0.74, that is over a distance comparable to 
one pressure scale height, which is about 
0.08, then the value of 70 would reach up 
to 0.8/xHz, that of 71 would be negative, 
reaching to -2.6, and that of 72 would also 
be negative reaching — 2.9/iHz. Clearly, 
such values are very significant and the 
field would be easily detectable. Detec- 
tion of the signal corresponding to a puta- 
tive 0.1 MG field is problematic at present 
day accuracy. The best chance is to see it 
is in the even-a coefficients, if indeed the 
field were dominated by the low-j poly- 
nomials. If the 1 MG field were present 
only within the overshoot zone, extending 
from (say) x=0.65 to the base of the con- 
vective zone, then the corresponding ex- 
treme values would be -0.31, -0.95, and -1.1 
//Hz. Thus, somewhat stronger than 0.1 
MG fields are required for detection. How- 
ever, stronger fields may be anticipated in 
the overshoot layer. 

Chou and Serebryanskyi (2002) found 
evidence for a 0.4-0.7 MG field at the 
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Fig. 7. — Kernels for calculating the val- 
ues of 7t,o according to eq. [73] , and aris- 
ing from dipolar (j = 1) and quadrupolar 
(j = 2) toroidal magnetic fields, for four 
£ = 30 modes of indicated orders n in the 
zone near the bottom of the convective en- 
velope {x = r/i?0 = 0.71). 
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Fig. 8. — The same as Fig. 7 but for the 
even-a splittings (A; > 0). 



base of the convective envelope from time- 
distance seismology. Such a field, if it per- 
sists over a distance to comparable to that 
assumed by us, should be detectable by 
means of global seismology. However, the 
effort made so far did not result in the de- 
tection of a significant signal (Basu, 2002). 
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Fig. 9. — Kernels for calculating the 7's 
arising from small-scale horizontal mag- 
netic field near the base of the solar con- 
vective zone. 

It is possible that the magnetic field in 
the deep part of the convection envelope, 
and in the overshoot zone forms azimuthal 
ropes, and thus is better represented 
small-scale field with its mean value being 
a slowly varying function of latitude. In 
this case, the frequency perturbation is de- 
scribed by the adiabatic version of eqs. [62], 
at /c = 0, and by [65], at A; > 0, with only 
the horizontal components included. For 
the 7's, we use an expression that is simi- 
lar to that given in eq. [68] 
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^"-/^^^"(t^)' '''' 

with X = t/Rq 

TT A mHz 9„rT 

= 1.48 X 10"^ ^ ^^M,fc- 

With the approximation for the eigenfunc- 
tions, which is valid for p-modcs in the 
zone considered as given in Eqs. [71] and 
[72] , we have 



Both kernels change sign in the region of 
interest. In Fig. 9, we show examples 
of the kernels for the same modes and in 
the same layer, as in Fig. 8. Differences 
between the figures are apparent. Note 
in particular, the sign changes within the 
layer. With a IMG field in the layer, we 
get 7o of about 0.8 /xHz for the n = 8 and 
10 modes. The absolute values at A; > 
are somewhat lower. So that a 0.4-0.7 MG 
field should be easily detectable in the fre- 
quency shifts. The expected signal in the 
even-a's would be somewhat weaker, but 
should also be detectable. 

10. Conclusions 

We surveyed various effects that may 
explain the increase of the mean frequen- 
cies of solar oscillations and the changes in 
the fine structure of the multiplets corre- 
lated with the rise of activity. Beyond any 
doubt is the fact that the seismic changes 
reflect temporal and longitudinal evolution 
of the sun's activity and that the main 



part of these changes has its source close to 
photosphere, perhaps within 1 to 1.5 Mm. 
The question that remains to be answered 
is how the changes are related to mag- 
netic field variations and associated varia- 
tion of the velocity field and the tempera- 
ture in the atmosphere, and in the subpho- 
tospheric layers. The answer is important 
because only after we know it, we will be 
able to make a full use of seismic data as 
a probe of the large-scale variability of the 
sun over its activity cycle. 

Another open question is whether cur- 
rent frequency data reveal changes in the 
deep interior. Assessing the intensity of the 
magnetic field in the lower convective en- 
velope and in the outer radiative interior, 
just beneath, would be of great importance 
for understanding the physics of magnetic 
activity. 

With these questions in mind, we devel- 
oped integral formulae expressing shifts in 
oscillation frequencies in terms of changes 
in magnetic and velocity fields and tem- 
perature. We made certain approxima- 
tions, which are well-justified both for so- 
lar p- and f-modes. Small scale-fields were 
represented by the square-averaged val- 
ues of vertical and horizontal components. 
These values were expressed in terms of the 
Legendre-polynomial series. We showed 
the coefficients at P2fe in these series con- 
tribute only to a2k coefficients of the fre- 
quency splitting. The connecting expres- 
sions were given in the form of integrals 
over the depth with explicit expressions for 
the kernels. We also provided kernels hnk- 
ing the mean frequency shifts and the a2k 
coefficients to a low-order expansion of a 
putative large-scale toroidal field near the 
base of the convective envelope. 
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Plots of the kernels allowed us to make 
a simple assessment of the changes needed 
to explain the measured shifts. We found 
that the increase of the mean frequencies 
and the changes in a2fe-coefficients are most 
easily explained in terms of a decrease of 
turbulent velocities associated with the in- 
crease of the magnetic field with grow- 
ing activity. The required decrease in the 
turbulent velocity needed to explain the 
data constitutes only a fraction of a per- 
cent. A decrease in turbulent velocity is ex- 
pected to result in a temperature decrease 
in outer layers of the sun. Our estimate 
showed that the resulting temperature de- 
crease should give a significant contribu- 
tion to the mean frequency increase, which 
reduces requirement on the size of the de- 
crease of the turbulent velocity. Account- 
ing for the seismic changes by the sole di- 
rect effect magnetic field is more difficult. 
An increase of the surface-averaged r.m.s. 
value of the vertical field component by 
about 0.1 kG between the minimum and 
maximum of the activity would be needed 
to account for the mean frequency increase 
of the p-modes. The measured changes in 
the even a coefficient for p-mode require 
about twice as large field increase at most 
active latitudes. Also a larger field seems 
to be needed to account for the systematic 
increase of the f-mode frequencies. 

Considering the influence of the fleld 
near the base of the convective envelope, 
we found that there is a chance for de- 
tecting such a field directly from the fre- 
quency data. Evidence, from time-distance 
seismology, for the presence there of the 
(0.4 — 0.7) MG field was recently put for- 
ward by Chou & Serebryanskyi (2002). 
We showed that such a field, if extends 



over a layer of thickness comparable with 

one pressure distance scale, should be de- 
tectable also by means of global seismol- 
ogy- 
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11. Kernels for the 7's 

Here we summarize the expressions for 
the kernels for evaluating the 7's due to 
small-scale velocity and magnetic fields. 
We begin by recalling the definition of the 
radial eigenf unctions: 

^ = r[yir)er+z{r)VH]Yrie, 4>) exp(-iu;t). 
We also define 

2 2 

. gr uj r 

and 

The kernels for centroid shifts (70) due 
to turbulent velocity are 

<isoth = Asoth - (A^ + 2A^A + A£) (Al) 

and those due to a small-scale magnetic 
fields are 

^M,isoth = 2A(zA + 2£) - Asoth (^3) 



^M,isoth = 2AVA(4^A+3£:)+Pi,oth. (^4) 
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where 



isoth 



-r 



2AzX 



Xp 

The kernels for the even-a coefficients 
(7fe, k > 0) due to turbulent velocity are 

7^];'fc = -(2C-V + A' + 2A^A + A£), (A5) 



T^v,k = C-2[^-r(i+rp)A2+2A^A+3A£:]. 

{A6) 

and those due small-scale magnetic fields 



where 



and 



4C-2V'-r(i+rp)A2+2A£:, (^7) 

^M,fe--2C + V' + 2A^ (^8) 



2 

Oi) r 



C= {yX + 2Azy) 



/J amp 



2 
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